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Abstract 

Inspired by theories such as Loop Quantum Gravity, a class of stochas- 
tic graph dynamics was studied in an attempt to gain a better understand- 
ing of discrete relational systems under the influence of local dynamics. 
Unlabeled graphs in a variety of initial configurations were evolved using 
local rules, similar to Pachner moves, until they reached a size of tens of 
thousands of vertices. The effect of using different combinations of local 
moves was studied and a clear relationship can be discerned between the 
proportions used and the properties of the evolved graphs. Interestingly, 
simulations suggest that a number of relevant properties possess asymp- 
totic stability with respect to the size of the evolved graphs. 

1 Introduction 

The dynamics of quantum-gravitational theories such as Loop Quantum Grav- 
ity can be formulated as local moves operating on a special class of labeled 
graphs. It would be interesting to know whether or not a similar set of local 
moves operating on generic graphs will generate an ensemble of graphs similar 
to that assumed by the such quantum-gravitational theories. To that end, unla- 
beled graphs in a variety of initial configurations were evolved using generalized 
Pachner moves. Theories of space-time generally assume underlying graphs 
which are, in an appropriate sense, approximately manifold-like. As it turns 
out, however, stochastic application of local moves does not seem to generate 
graphs which are approximately manifold-like. Although the graphs produced 
do not seem directly interpretable as manifold-like structures, simulations sug- 
gest that the ensemble of such graphs is distinguishable from the ensemble of 
random graphs with the same size and average degree. Certain properties of the 
graphs produced by the evolution clearly depend on the proportions of different 
local moves used to generate the evolution. Additionally, simulations suggest 
that these same properties possess asymptotic stability with respect to the size 
of the evolved graphs. 
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2 Definitions 



A graph is defined to be a pair G = (V, E), where V is a set of vertices, and 
E is a set of vertex pairs called edges. The degree of a vertex is the number of 
edges incident on that vertex. The graphs described in this paper are unlabeled 
and contain no self-loops (edges connecting a vertex to itself) or duplicate edges 
(edges connecting the same two vertices as some other edge). 

2.1 Generalized Pachner Moves 

For planar trivalent graphs, some combination of the three Pachner Moves can 
convert any graph into any other graph. The graph dynamics discussed here is 
a combination of three distinct transformations which are each generalizations 
of one of the three Pachner Moves. As with traditional Pachner moves, one 
is an expansion operation, one a contraction operation, and one an exchange 
operation. The expansion move involves replacing a vertex with a d-simplex (A 
d-simplex is a set of d vertices where each vertex is connected to all of the other 
d—1 vertices, i.e. a triangle for d = 3) (Figure QJ. The contraction move is the 
inverse of the expansion move (Figure [2J ■ The exchange move exchanges the 
endpoints of edges both adjacent to a third edge (Figure |3J). Specifically, given 
some edge E\ = (Vi, V 2 ) and edges E 2 = (Vi, V3) and E 3 = (V2, V4) replace E 2 
and E3 with E' 2 = (1/1,14) and E' 3 = (X^, V3). Note that the exchange moves 
do not preserve planarity or any manifold- like structure (Figure 0J. 



Figure 1: Expansion ("1-to-d") Pachner Move, d = 3 is shown here. 



When applying an expansion move, vertices with a degree less than d were 
not considered for replacement. The old incoming edges were distributed among 
the new vertices such that each new vertex in the inserted d-simplex was assigned 
at least one of the old edges. 

2.2 Dimension 

The statistical dimension of the produced graphs was computed. Let d(v, w) be 
the number of edges in the shortest path between vertex v and vertex w. Let 
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Figure 2: Contraction ("d-to-1") Pachner Move, d = 3 is shown here. 




V(v, r) be the number of vertices w such that d(v, w) < r. If, for a generic vertex 
v, V(v, r) = Cr a for some a then a is the statistical dimension. It should be 
noted that, for an arbitrary graph, the statistical dimension, as specified here, 
is not well defined. 

2.2.1 Analytic Computation 

It is possible to compute the function V(v, r) for some generic vertex v in a 
connected random graph with n vertices and m edges. 
Let d be the average degree of a vertex. Then: 

d=^ (1) 
n 

The vertices in the graph can be divided into equivalence classes based on 
the length of the shortest path to vertex v. For any given path length r the size 
of the equivalence class will be denoted A(r). It follows immediately that: 

V(r) = V(r-l) + A(r) (2) 
and also that the number of vertices not within a distance r — 1 of v is: 

R(r) = n-V(r-l) (3) 
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Figure 4: Another Exchange Move 




s(r) is the average number of "free edges" per vertex at a distance of r. 

s(r)=f(r-l)/A(r-l) (4) 
A{r) is the number of vertices at a distance of r. 

A(r)=R(r)(l-(l-^) A ^) (5) 

f(r) indicates the number of "free edges", meaning the number of edges 
connecting vertices at a distance of r to vertices at a distance of r + 1. 

f(r) = (dA(r) f(r 1))(1 - 1(^-^1)) (6) 

Plots of this relation for various values of n and d are shown in Figure |S] It 
is interesting to note that a generic feature of the random connected graphs is a 
small statistical dimension and small distances and a larger dimension at larger 
distances. 

2.3 Non-local Edges 

An edge is said to be "non-local" if it is not part of any short cycle (where short 
means the length of the cycle is less than or equal to some rj). Intuitively, this 
should mean the edge connects locations in the graph which are "far away" . If 
the number of non-local edges is small compared to the total number of edges, 
then it may be possible to view the graph as some "local" structure (such as a 
regular lattice) combined with some additional distribution of edges. 

3 Description 

3.1 Input Parameters 

The simulations used to investigate the properties of the chosen graph dynamics 
depended on only a few input parameters. These parameters were: 

• d: The number of vertices in the simplex added during an expansion move. 



4 




• n : The length of the largest cycle to consider local. 

• R: The ratio of exchange moves to expansion moves. 

• R c : The ratio of contraction moves to expansion moves. 

3.2 Graph Pictures Coloring 

The graph pictures which appear in this paper were drawn using an algorithm 
which treats each edge as a spring with some spring constant and solves for the 
equilibrium configuration. 

The edges are colored as follows (where r; = 5): 

• Red - Initial edges which are "local" 

• Magenta - Initial edges which are "non-local" 

• Black - Edges added by expansion moves which are "local" 

• Green - Edges added by expansion moves which are "non-local" 
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3.3 Initial Conditions 



Many different starting configurations were tested including a single triangle, 
regular lattices, almost complete graphs ( d = (1 — e) * (n — 1)) and random 
connected graphs. The exact starting configuration seemed to have very little 
long-term influence on the simulation results. An example of a random con- 
nected graph used as an initial configuration is pictured in Figure 



Figure 6: Random Connected Graph (n — 200 and m = 400) 




3.4 Simulation Outputs 

For the purpose of this study, the graphs were evolved until they contained tens 
of thousands of vertices (except for cases when R c = 1). Unfortunately, due to 
the slow convergence of the graph-drawing algorithm, graphs of that size could 
not be rendered in a practical manner. However, the properties of the graphs 
being evolved often converged well before reaching several thousand vertices, 
and graphs of that size could be feasibly rendered. Figure [7] shows the result 
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of the graph pictured in Figure evolved using R — to where n = 2200 and 
m = 3400. Figure El shows a portion of the dense center region and Figure 
shows the end of one of the spikes. Approximately 14.7% of the edges are "non- 
local." Figure EH shows the result of the graph pictured in Figure evolved 
using R = 100 to where n = 1800 and m — 2800. Figure ITT1 shows a portion 
of a central region. Approximately 42.5% of the edges are "non-local." Notice 
that the addition of the exchange moves seems to suppress the formation of the 
spiky structures. 

Figure 7: Evolved Random Graph (Evolved using R = to n = 2200 and 
m = 3400) 
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Figure 8: Center Region of Evolved Random Graph (Evolved using R = to 
n = 2200 and m = 3400) 




4 Results 
4.1 Dimension 

The graphs which resulted from the evolution bore little resemblance to regular 
lattices or any other structure with a statistical dimension constant with scale. 
V(r) for these graphs is qualitatively similar to that for a random connected 
graph as can be seen in Figure Upland Figure 1131 which show \n(V(r)) vs. ln(r) 
for an evolving graph with R = 1 and R = 100 respectively (a regular lattice of 
dimension d would appear as a line of slope d in those figures). 

Noting that the dimension on small scales is closer to 1 in Figure IT51 where 
R = 1 compared to FigurefHflwhere R = 100. This can be understood by noting 
that the graphs evolved with R = 1 were dominated by spiky structures which 
are nearly one-dimensional, whereas evolution with larger values of R suppressed 
the spike formation. FigurelTHshows the approximate slope of ln(V(r)) for small 
r (as determined by a heuristic algorithm) vs. R. After some value of R the 
dynamics appears to become "exchange dominated," whereby further increasing 
R does not significantly change the structure of the graph. This can also be 
seen using other measures (such as graph diameter). 
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Figure 9: Spike of Evolved Random Graph (Evolved using R = to n = 2200 
and 77i = 3400) 



4.2 Non-Local Edges 

Statistical dimension and similar measures do not completely characterize the 
structure of a graph. The percentage of edges which are "non-local" can provide 
an indication of how much the graph resembles a "local" structure such as a 
regular lattice. Figure ^| shows a plot of the fraction of non-local edges vs. 
n for some graph evolution (for r; = 5). For any value of R it is possible to 
identify the asymptotic value of the fraction of non-local edges. FigureHBlshows 
the fraction of non-local edges vs. r\ for some fixed graph size. It should be 
noted that the fraction has only reached its asymptotic value for n less than 
the graph diameter. Figure IT71 shows the fraction of non-local edges vs. n for a 
graph evolving using R = 100. 

5 Conclusion 

The properties of unlabeled graph evolution under the action of local moves have 
been investigated and several interesting properties seem to be asymptotically 
stable with respect to the number of vertices of the graph. Evolving the graphs 
with < R c < 1 did not significantly affect the long-term results (larger values 
of R c , as should be expected, result in graphs of static or decreasing size). 
Also, evolving the graphs with d > 3 results in qualitatively similar behavior 
(although, for example, the asymptotic fraction of non-local edges changes as a 
function of d) . 
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Figure 10: Evolved Random Graph (Evolved using R = 100 to n = 1800 and 
m = 2800) 
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Figure 13: \n(V(r)) vs. ln(r) for graph evolved with R = 100 
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Figure 14: Small-scale dimension vs. R 
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Figure 15: The fraction of non-local edges vs. n (77 = 5) 
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Figure 16: The fraction of non-local edges vs. r; (n = 5400) 
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